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We study the distillability of a certain class of bipartite 
density operators which can be obtained via depolarization 
starting from an arbitrary one. Our results suggest that non- 
positivity of the partial transpose of a density operator is not 
a sufficient condition for distillability, when the dimension of 
both subsystems is higher than two. 
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I. INTRODUCTION 

Maximally entangled states represent an essential in- 
gredient in most applications of Quantum Information 
(QI) |l|]. In particular, in Quantum Communication one 
can use them for transmitting secret messages between 
two locally separated parties [g[. In practice, however, 
states are mixed due to the interaction with the envi- 
ronment, and are not usable for those applications, even 
though they may be entangled. The solution to this prob- 
lem was presented by Bennett et a/., Deutsch et al., and 
Gisin H-ra] , who have given a procedure to "distill" max- 
imally entangled states of two qubits out of a set of pairs 
in certain (mixed) entangled states, by only using local 
actions and classical communication Mm- Later on, the 
Horodecki family showed that any, even infinitesimally 
entangled state of two qubits (two-level systems) can be 
distilled into a singlet B . They have also proved a neces- 
sary condition for the state of an arbitrary bipartite sys- 
tem to be distillable, namely, that the partial transpose of 
the corresponding density operator must be non-positive 

As shown by Peres |10|, the positivity of the partial 
transpose is a necessary condition for separability. In 
fact, this condition turns out to be a sufficient condi- 
tion for separability in both cases of qubits (two-level 
systems), or one qubit and one trit (three-level system) 
Ipdl . A natural question arises: is this condition also suf- 
ficient for separability for higher dimensional systems? 
P. Horodecki has recently shown |g2| that there are in 
fact states in higher dimensional systems which have a 
positive partial transpose, but are non separable (see 
also P,|l3|-[l7| ) . As a consequence, positivity of a par- 
tial transpose is, except for 2x2 and 2x3 systems, 
not sufficient for separability. Similarly, non-positivity of 
a partial transpose is necessary for distillability, and is 
sufficient for 2x2 and 2x3 systems. However, the ques- 



tion whether this condition is sufficient for distillability 
in higher dimensional systems remains still open. 

In this paper we investigate distillability of high dimen- 
sional systems shared by two parties, Alice and Bob. We 
introduce a depolarizing superoperator that allows one 
to reduce an arbitrary density operator with non posi- 
tive partial transposition (NPPT) to one with the same 
property, but in a standard form that is characterized by 
a single parameter. We analyze some properties of those 
operators, and show that for any given finite number of 
copies there are density operators p for which one can 
never find a subspace of dimension 2 in Alice and Bob's 
Hilbert spaces in which p still has NPPT. We also present 
some numerical evidence that indicates that this class of 
states is independent of the number of copies. All these 
results suggest that there exist states with a NPPT and 
which are not distillable. 

This work is organized as follows: In Section |l| we 
review some of the present knowledge concerning distil- 
lability and entanglement, and we introduce the defini- 
tions and properties that are needed in order to study 
the problem of distillability of general density operators. 
In Section HI we concentrate on the case in which Alice 



and Bob have three-level systems, whereas in Section EV 
we generalize our results to the d-level system case. In 
Section M we show the basics of the numerical procedure 
used to study the distillability of 2 and 3 copies. Finally, 
we summarize our results. 



II. ENTANGLEMENT AND DISTILLABILITY 

We consider two parties, Alice and Bob, who share sev- 
eral pairs of particles. Each pair is in a state described by 
the same density operator p. We will assume that Alice's 
(Bob's) particles are d^-level systems (ds); that is, the 
density operator p acts on the Hilbert space (p ^ ®(p ^ . 
We will denote by {|1), |2), . . . , \dA)} an orthonormal ba- 
sis vii(p ^ and analogously for (p ^ . We will also use the 
notation \i,j) = \i) a ® |i)s- 

We will assume that Alice and Bob are able to manipu- 
late their particles by only using local actions (operators 
and measurements) and classical communication. In this 
case, we say that the density operator p is distillable if 
they can produce a maximally entangled state 



l<I>rf) = 



1 '^ 



Vd 



(1) 



where d = niin(d^, ds)- On the other hand, we say that 
p is separable if it can be prepared out of a product state 

(e.g. |1,1)). 

In this Section we will review some of the results de- 
rived by Peres and the Horodecki family concerning dis- 
tillability and entanglement, and will introduce the defi- 
nitions and properties that are needed in order to study 
the problem of distillability of general density operators. 



A. Partial transposition 

As shown by Peres 0] and the Horodeckis |§,|ll|,|l2|, 
the partial transpose of a density operator plays an im- 
portant role in establishing its distillability and entan- 
glement properties. In general, given an operator X 
acting on ^ "* (g) ^ ^ , we define the partial transpose 
of X with respect to the first subsystem in the basis 
{|l),|2),...,|dA)},X^^, as follows: 



XTa 



i,j = l k,l=l 



't,k\X\jJ)\j,k){tJ\ 



(2) 



In the following we will use a property of this operation, 
namely tr{YX^^) = tr(y^'^X). 

We say that a self-adjoint operator X has a non- 
positive partial transpose (NPPT) if X^-^ is not positive; 
that is, if there exist some |*) such that {^\X'^'^\'i') < 0. 
The positivity of the operator p^-* gives necessary cri- 
teria for separability and non-distillability of a density 
operator p. In particular: (1) If p is separable, then 
P^^ > @; (2) If p^-* > then p is not distillable 
|l2[ . These two necessary conditions turn out to be suf- 
ficient for dA = 2 and ds < 3 g]. However, it has been 
shown that the first condition is not sufficient for separa- 
bility for the rest of the cases {dA — 2 and ds > 3, and 
dAjds > 2) |12|. On the other hand, nothing is known 
about whether the second condition is also sufficient for 
non-distillability in these cases. 



B. Distillability 

The problem of distillability of general density opera- 
tors acting on ^ "^ ®(p ^ can be expressed in a simpler 
form |18] . A density operator is distillable iff for certain 
positive integer N, we can find a state of the form 



such that 



\'i)=a\e,)A\fi)B+b\e2)A\f2)B, 



(*|(p®^)^^|*)<0. 



(3) 



(4) 



Here, {|ei)y!i, |e2)yi} are two orthonormal vectors in 
(^ '^)®^, and {|/i)_B, 1/2)3} are two orthonormal vec- 
tors in {(p ^)®^. This condition basically means that if 



Alice and Bob share N pairs, one just has to find a two- 
dimensional subspace in the whole Hilbert space of Alice, 
and another in Bob's such that the projection of p®^ in 
such subspaces has NPPT. The reason is that if one finds 
such a subspace, then according to what was exposed in 
the previous subsection one can distill a maximally en- 
tangled state in 0" ^(p , which can be converted into 
a maximally entangled state in ^ ^ ®(p ^ . Conversely, 
if one can create one of those states then one can also 
produce one in ^ (E)(p , and therefore this ensures that 
(|) must be fulfilled. 

Thus, in practice, one can analyze for each number of 
copies A^ = 1, 2, . . . whether the condition (Q) is fulfilled. 
In order to facilitate this task, we will use the following 
definitions: If for a given N condition (Eh is fulfilled we 
will say that p is A^-distillable. On the other hand, if 
for a certain N there does not exist any Ivf) satisfying 
Eq. (Eh, then we will say that p is A^-undistillable. Thus, 
p is distillable iff there exists an A^ for which it is A^- 
distillable. Conversely, p is non-distillable iff it is A^- 
undistillable \/N. 



C. Distillability in(?^ (g,(/;'^ 

With the properties and definitions given above, one 
can very easily prove that when dA — '2 and ds > 2 and 
if p has a NPPT then p is 1-distillable. The reason is 
that there exists some JVP) such that (»I'|p'^-*|\I') < 0. On 
the other hand, since d^ = 2 then the Schmidt decompo- 
sition of 1^) has at most two terms, and therefore can be 
written in the form (0). Thus, in this case non-positive 
partial transpose of p is a necessary and sufficient condi- 
tion for distillability. 



D. Depolarization in^ ®(? 

In this subsection wc will introduce some superoper- 
ators which will be useful to study if a given density 
operator is A^-distillable. We will also show that given 
a density operator one can reduce it to a standard form 
which is characterized by a single parameter [H, and 
that preserves the distillability properties of the original 
state. 

Let us first define some useful projector operators. 
Given two quantum systems with corresponding Hilbert 
spaces (p , we denote by n^ the permutation operator, 
and by 



Ad = il- Hd)/2, Sd = t-A, 



nd)/2, (5) 



the projector operators onto the antisymmetric and sym- 
metric subspaces of ^ (p , respectively. Note that 
tr(S'd) = d{d + l)/2 and tr(Ad) = d{d ~ l)/2. We also 
define the projector operators 



p. 



^d){M. Qd = i-Pd 



(6) 



where |(/)d) is the maximahy entangled state defined in 
(0). One can easily check that 

Pj- = i(l-2A,), Al- = \{l-dPa). (7) 

We define the depolarization superoperator T), acting 
on any given operator X, as follows: 



^W = A.^I(#^ + 5.M^^^) 



tr(Ad) 



tr(5d 



(8) 



This superoperator is a projector, is self-adjoint (on the 
Hilbert-Schmidt space of operators acting on 0" ®(p ), 
and preserves the trace. In the Appendix A we show that 
we can write ilSl 



V{X)^ dnu{U(g>U)X{U(g>U)'' 



(9) 



where the integral is extended to all unitary operators 
acting on ^ and / d/iu = 1 [dfJ- represents the standard 
invariant Haar measure on the group SU{d)]. We will 
later on use the partial transpose of T>{X), and to this 
aim we define the following superoperator 

£{X) = [V{X^^)f^ = / d^u{U* ® U)X{U* ® U)^ (10a) 

triQdX) 



PMPdX) + Qd- 



tr(Qd) 



(10b) 



where U* denotes complex conjugation in the basis in 
which the partial transposition is defined. This superop- 
erator is also a projector, self-adjoint, and preserves the 
trace. Note that for any unitary operator V acting on (p 
we have 



{V ® V)V{X){V ® V)^ = V{X), 
{V* (g) V)£{X){V* (g) Vy = £{X). 



(11a) 
(lib) 



The form (^ shows that the superoperators V and £ 
can be implemented by means of local operations. In 
particular, it shows that any density operator p can be 
transformed, using local operations, to the form 



T^ip) = Pa 



Nia 



■{Sd + aAd) 



(12) 



where a is such that ti^Adp) — tT{AdPa) and N{a) = 
tr(S'rf) + atT(Ad) is a normalization constant. That is, 
one can depolarize any density operator to the one pa- 
rameter family (|l2| ) while keeping the weight in the an- 
tisymmetric subspace. We will be more interested in the 
partial transpose of pa , which is given by 



[V{p)Y-=£{p'-) 



Ta-^ 



1 



M(/3) 



[Qd'PPd), 



(13) 



where M{(3) — tr{Qd) — /3 is a normalization constant, 
and the relationship between a and /3 is /3 = [{a— l){d — 
1) - 2]/(a + 1) < d - 1. Note that since 



Pa=2?(Pa), pi^^Sipi^) 



.Ta^ 



we have that 



{U* ®U)pI^{U* (^U)^ ^ pI^. 



(14) 



(15a) 
(15b) 



for any unitary operator U . 

Using the properties derived above, one can easily 
check that 



Pa is separable <^ p^'* > <^ /? < 



(16) 



The last equivalence follows directly fro m ([l3|) . For the 
first one we have: (=>) See subsection ( [I A ); (<;=) We 
have that for (3 — [i.e. a = a^ = {d + l)/(d — 1) = 
iv{Sd)/iy[Ad)] 



p^;^aQ<iCx£(|0,l)(0,l|), 



(17) 



which is obviously positive and separable, in which case 
the same holds for pag- For /3 < (a < ao) we can 
always obtain pa by adding the identity operator (which 
is separable) to pag (this is due to the fact that tr(S'c() > 
iy[Ad)). 

Thus, for /? > (a > ao) Pa is non-separable. One 
can easily check that this condition is equivalent to 



tr{AdPa) > 1/2. 



(18) 



This last form allows us to show that for any given 
density operator p with NPPT, one can always trans- 
form it using local actions to the form (12) such that 
it still has NPPT ^. Let us show that. Suppose 
that for a given |*), (*|p^-*|vl') < 0. We can write 
I*) = Z]r=i Ci\u,,Vi) where \ui)^^-^ and \vi)^^-^ form an 
orthonormal basis. The operator p can be transformed 
by local operations to ps ex {A'' (g) B'')p{A (g) B), with 
tr(ppP„) < (by simply taking A = Y:7=i l<)(*l/c. + 
Etn+il<)(*l and B = Et„+il^''^>(*l)- Using that 
> tr(p^''P„) = tr(psPj-'i) and (^ we immediately ob- 
tain that tr{AdPs) > tr(yl„ps) > 1/2. Since V conserves 
this quantity, we obtain that 'D{ps) has a negative partial 
transpose. 

As pointed out by the Horodecki iQ, the problem of 
distillability can be reduced to the study of density oper- 
ators of the form (12). If we find that all those operators 
with a NPPT are distillable, we will have shown that 
NPPT is a necessary and sufficient condition for distil- 
lability. On the contrary, if we find that there exist an 
operator of the form ([l^ ) which has a NPPT, but is not 
distillable, we will have shown that such a condition is 
not sufficient. In the following Sections we will show that 



there are density operators of the standard form with a 
NPPT which are not 7V-distillable for certain values of 
N. As we have seen in subsection II B we can study that 



by checking whether there exist vectors of the form 
fulfiUing condition (m. 



III. iV-DISTILLABILITY IN (?^ (g)(?^ 

We consider the case d^i = d^ = 3 and a density oper- 
ator of the form 



1 



N{a) 



(S + aA), 



(19) 



where we have omitted the superindices d = 3, and 
N{a) = 6 + 3a. According to the discussion in Sub- 



section (IID) we just have to consider a > uq = 2, since 
otherwise pa is separable. We also have for the partial 
transpose 



P/3 



1 



M(/3) 



(0 - f3P) 



(20) 



where /3 = (2Q!-4)/(a + l), with 2 > /3 > and M(/3) 
8-/3. 



A. 1— distillability 

We look for a vector of the form (||) such that ^ is 
fulfilled. Cho osing U such that U\ei^2) = |1, 2) and using 
the property (|l5b|) we see that we can restrict ourselves 
to the subspace spanned by {|1), |2)} in Alice's Hilbert 
space. Defining by I2 the projection operator into this 
subspace and by I3 the identity operator in (p , we obtain 
after projecting p^'^ onto such subspace 



l^pfl^^tt 



2(1 +/3) 



P2 



(21) 



which is positive iff /3 < 1/2. Thus, we obtain that Pa is 
1-distillable iff /? > 1/2 (or, equivalently, a > 3). 



B. 2-distillability 

Let us consider now two pairs, in a state pa- We will 
show that for (3 < 1/4 the state pa is 2-undistillable. For 
any state | ^) of the form (0) we have 



tr2[tri(|vI/)(M/|g(i))(Q(2)-p(2)/2)]. 



(22) 



where the superscripts 1 and 2 refer to the first and sec- 
ond pair, respectively. Using the fact that Q is sep- 
arable, and therefore that it can be written as Q = 
X]jCi|ai,6i)(ai,5i| with a > 0, we have 



tri(|vl')(vI/|QW) = ^c.|*,)(*. 



(23) 



where \'i'i) — {ai,bi\'^) is a state acting on the sec- 
ond pair which itself has the form (0). Thus, accord- 
ing to the results of the previous subsection we have 
that ('J'i|(Q(2)_p(2)/2)|^^) > and therefore (*|Q(i)(^ 
(Q(2) - P(2)/2)|$) > 0. In the same way we have that 



0<(*|(Q 



(1) 



p(i) 



)®(Q 



(2) 



p(2) 



) - P(l) 



p(2) 



m 



(24) 



Using the fact that P*^^-' (g) P^^-* > we obtain the desired 
result. Note that our results do not imply that for 1/4 < 
(3 < 1/2 Pa is 2-distillable. In fact, as shown in the 
next section, numerical calculations indicate that it is 
2-undistillable. 



C. 7V-distillability 

We consider now N pairs, in a state pa- We will show 
that for f3 < 4~^ the state pa is iV-undistillable. For any 
state IvP) of the form (p|) one can check the following rela- 



tions: (i) (^IP®'^!*) < ^, (ii) (*|Q 



and (iii) (*|Q»^|*) 



> 



3fc J 

1 



^N-k 



p®/c|^^ < Jr 



To show (i), one uses P^ 



P(ik and the property (15b), from which follows that the 
projection into the subspace spanned by {|1), |2)} gives 
the maximum value for {'^\P'^^Y^) . From (i) we immedi- 
ately obtain (ii) by using that (*|QX|*) < {^tX\^) for 
all positive operators X . Relation (iii) can be obtained 
by using (i) and the separability of Q in a similar way as 
in the previous Section. Combining (ii) and (iii) we find 



makQ 



®N 



Eq' 



-k TjCSk 



p««|^^ >0, 



(25) 



with afc = 2 ( ^ ) 3^ ^ and the sum runs over all possible 
permutations of the pairs. By summing (P5|) for all odd 



k and using that X^feodd 
for /3jv < (3n 



Qk ^ 4: 



N 



■)N - 



Pff , one finds 



< (*|g®^ -PnY^Y. Q®^"''P®''|*) 

k odd perm 
fe odd porm 



<(vI/|(Q-/3jvP)®^|*). 



(26) 



We used that maxfe(/3^) < Pn (for (3n < Pn < I) to 
obtain the first inequality (line 2), while we added all 
positive terms (even k) in the second step (line 3). This is 
already the desired bound, i.e. for < /3 < 4~^ < /^^r = 
jT^^^w, the state pa is A^-undistillable. Again, this does 
not mean that pa is A^-distillable for 4~^ < /? < 1/2. 
In the Appendix C we present a better bound for (3. 



IV. iV-DISTILLABILITY IN (p° 



A= (*|i?|*) <0, 



(28) 



We consider the case (Ia = ds = d and a density oper- 
ator of the form ( |l2| ) with the partial transposition given 
by®. 

Similar techniques as in the d — S case can be used to 
obtain bounds also for arbitrary d. One finds for example 
that Pa is 1-distillable iff /3 > | — 1. We also obtain that 
Pa is 2-undistillable for /? < ^^ and iV-undistillable for 

/? < min(/3Ar, /3^'^) with Pn = (d+iy^ {d~i)^ ■ ^°*^ *^^* 
the minimum is required here, since - differently to d—3 
case - one can have that /3jv > 1. In this case one has to 
chose Pn < /3^/^ which implies maxfc/3^ < (3^ to ensure 
that the first inequality in (EQ) remains valid. 

Furthermore, there is an interesting relation between 
the states pa for different d. Imagine we would like to 
convert a single copy of a state pa in ^ (g)^ to a state in 
^ ®(p {k < d) in an optimal way, i.e. to obtain a new a 
which is as large as possible. We show here that whenever 
we convert a state pa to some lower dimension, there will 
always be some states which loose the negativity of their 
partial transposition. In order to prove this, we consider 
vectors |\I'fe) with k Schmidt coefficients and show that 
{^d\pl^\^d) < 0, while (*fc|pl^|*fc) > OV|*fc). Due to 
the property ( |l5bD , one can restrict oneself to the sub- 
space spanned by {|1) . . . |fc)}. Let us denote the identity 
operator in this subspace by 1^. One finds after project- 
ing p^-* onto Ife in A and B 

h «> Ifc Pa h ® Ifc « Ifc ® Ife ^ Pk (27) 

which is positive iff /3 < f — 1, while p^-* before the 
projection was positive iff /3 < 0. Thus all states with 
< /9 < f ~ 1 los'5 the negativity of their partial transpo- 
sition after the optimal projection onto a fc-dimensional 
subspace. The new (3k can be calculated from f3d of the 
initial state by Pk = ^iPd + 1) — 1- 

Finally, let us consider N copies of pa of dimension d, 

which can be viewed as a state in ^ ® (p . With 
tr {AdPa) = Ad, one finds that the state 'D{p®'^) = pa - 
the state in the high dimensional Hilbert space after de- 

polarization - has X^n = \.Y{AdNpa) — — ^^ — ^ '' • O^*^ 
checks that for Ad > 1/2 (i.e. pa is inseparable) we have 
ArfN < AdVTV, which simply means that the weight in the 
antisymmetric subspace decreases when going to more 
copies. Note that using this notation, we have that the 
state Pa is separable for A < 1/2, while it is 1-distillable 
for A > 2(2d'-i) ' '^hich tends to | for d ^ oo. 



V. NUMERICAL PROCEDURES 

In general, we are interested in showing there exists a 
JVP) of the form (||) for which condition (^) is fulfilled, i.e. 



where we have defined R ^ {Q - /3P)®^ with < (3 < 
1/2. In order to check that, we can minimize A with re- 
spect to |ei^2), 1/1.2) and a while keeping the normaliza- 
tion and orthogonality relations. One can readily check 
that the minimization implies 



(eili?!*) 
(e2|i?|*) 

ifiim) 

(/2|i?l*) 



Aoa|/i), 
A0&I/2), 
Aoa|ei), 
Ao6|e2). 



(29a) 
(29b) 
(29c) 
(29d) 



Note that the operator A{ei\R\ei)A is strictly positive, 
and therefore invertible. The reason for that is that for 
any |/)b, we have (e, /|i?|e, /) = (e*, /|i?^^|e*, /) > 
since according to ( [1^ ) we can always write R^^ = cl+B 
(where B > and c > 0). On the other hand, a,5 ^ 
since otherwise |\1/) would be a product vector and there- 
fore A > 0. Thus, we can use (29a) to write 



1/1) = - 



1 



a Ao - (ei|i?|ei) 



(ei|i?|e2) 1/2) 



which, after substituting in Eq. (29b) gives 

(e2|F(Ao)|e2)|/2>=Ao|/2), 

where we have defined 

1 



F(Ao) = i? - i?|ei) 



(ei|i?|ei) 
The normalization of |/i) gives 

1 



Ao 



(eilR. 



= (e2,/2| 



(ei|i?|ei)-Ao 



1^2, /2)- 



(30) 



(31) 



(32) 



(33) 



Thus, the problem is reduced to showing whether Eq. 
( pi| ) possesses solutions for Aq < 0. In that case, we can 
find a, 6, 1/1) using the other equations. On the other 
hand, if we have that 



(e2|F(0)|e2)>0 



(34) 



for all |e2), then we will have that there exists no solution 
with Ao < 0. This is so since i^(-|Ao|) - F(0) > 0. 

We have made a systematic search of the states 
|ei,2) which minimize the minimum eigenvalue of 
(e2|F(0)|e2) > for d = 3 and A^ = 2,3. Note that 
for N — 2 copies we can simplify further the numerical 
search by using the symmetries of the problem, which 
imply that we can choose |ei) = X)i=i "^jKi *) with c; > 0. 

In both cases we have found that this minimum eigen- 
value is > for (3 < 1/2, which strongly indicates that 
Pa is 3~undistillable (and hence also 2-undistillable) for 
(3 < 1/2. This is exactly the same bound that we had 
obtained analytically for 1-undistillability. 



VI. CONCLUSIONS 

We have shown that in order to study the distillabihty 
properties of bipartite d-level systems, it is sufficient to 
consider only the one-parameter class of states pa (12). 
By investigating the distillability properties of this fam- 
ily of states, we found strong indications that this fam- 
ily provides examples for non-distillable states with non- 
positive partial transposition. In particular, we found 
that for any given number of copies N there exist A'^- 
undistillable states which have NPPT. Guided by the re- 
sults of the numerical investigations, we conjecture that 
for d = 3 the states pa are non-distillable for /3 < 1/2, 
while they have NPPT for /? > (see also Fig.|l|). 



Note added: After completing this work we became aware 
of the results of D. P. DiVincenzo et.al ||22l, in which they 
also found evidences for the existence of non-distillable 
states with NPPT. 
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APPENDIX A: INTEGRAL REPRESENTATION 

OF V{X) 

We show here that the superoperator V as defined in 
(0) can also be written in the form (0). We restrict our- 
selves here to operators X which are density operators p 
for convenience, but exactly the same line of arguments 
holds for arbitrary self adjoint-operators X. As shown in 
Appendix B, we have that the depolarization superopera- 
tor T) can be implemented by a finite sequence of bi-local 
operations (^) of the form U ®U . Furthermore we have 
that the projector onto the antisymmetric subspace is in- 
variant under unitary operations of the form U ®U , i.e 
U ®U AdW ®W — Ad, which can be easily seen by using 
that Ad — {I — Hd)/2. From this property automatically 
follows that also pa dl2) is invariant under unitary op- 
erations of the form U ®U , since po, ex (1 -I- aAd). It is 
now straightforward to show (0) by using for any V 



dnu{U(g)U)p{U(g>U)'' = 

dpu'iU' ® U'){V (g) V)p{V (g> V)''{U' (g> U')K (35) 



with U'V — U. Taking p^ such that J^kPk = 1 we can 
write 



dpu{U®U)p{U®U)'^ = 
J2pk f dpu'iU' U'){Uk ® Uk)piUk ^ Uk)\U' ® u'Y 



dpu'iU' (E) U')p„{U' ® U')'^ =Pc.= T){p), 



(36) 



where we used ( pS]) in the first equality, while the sec- 
ond equality follows from (^ and we finally used the 
invariance of pa under operations of the form U ® U . 
This already shows that (Q) is fulfilled, i.e. 2?(p) = 

Sdpu{U®U)p{U(gUy. 



APPENDIX B: DEPOLARIZATION 

We are going to show now that an arbitrary state p 
can be depolarized to the standard form (II^) by a finite 
sequence of bi-local (random) operations without chang- 
ing the weight in the antisymmetric subspace Ha{d), i.e 
we show that there exist unitary operators Uk and prob- 
abilities pk such that 



^PkUk<E) UkpUl (E) Ul = V{p) = Pa 



(37) 



with tr (Adp) ~ tr (AdPa)- We first introduce a basis of 



1^^) 



— {\^)A\J)B±\J)A\^)B) 



\Xk) - \k)A\k)i 



(38) 



with i < j and {i,j,k) G {!,..., d}. Note that |(/3^) pro- 
vides a basis of the antisymmetric subspace Ha{d), while 
{|(^^), IXfc)} is a basis of the symmetric subspace Tls{d). 
The projectors into the symmetric/antisymmetric sub- 
space can thus be written as 






Xl'^'^^^^'^l 



(39) 



fe=i 



Let us write p in the basis (pq). In order to prove the 
statement (p7\), we devide the depolarization procedure 
into three steps: 

• we show that one can make p diagonal in the basis 
(pq) without changing the diagonal elements. 



• we prove that the antisymmetric subspace Ti.a{d) 
can be mixed up, i.e. one can equalize the coeffi- 
cients of |</3j^)(</3Z I without changing the weight in 
Haid). 



• finally we show that also the symmetric subspace 
can be completely mixed up without changing the 
weight in Tis{d). 

These three steps together ensure that p can be depolar- 
ized to the standard form pa. 



A. Diagonalizing p 

By mixing we understand in the following that a cer- 
tain operation U is (randomly) performed with probabil- 
ity p by Alice and Bob, while with probability 1 — p no 
operation is performed. The resulting density operator 
after this mixing-operation reads 



Pncw = piU «) [/pC/t U^) + (1 - p)p. 
We define the operation Ui as follows: 



Ui\k) 



ciTrSki 



\k), 



(40) 



(41) 



i.e. the state |/) picks up a minus sign while all oth- 
ers remain unchanged. Let us perform a sequence of d 
mixing operations, using p~^ and Ui with I = {1, ...,d} 
respectively. One can easiliy check that all diagonal ele- 
ments remain unchanged, while all off-diagonal elements 
of the form \ipfj){xk\ and \ipfj){^fij,\ for {i ^ i' , j ^ f) 
are eliminated. 

Let us define now the operation Ui , which introduces a 
phase i for the state |^) while it leaves all other states un- 
changed. Performing again the same sequence of d mix- 
ing operations as described above, but with Ui replaced 
by Ui , one can check that all off-diagonal elements of the 
form |xi)(Xjl are eliminated. 

We finally define the operation Uki, which simply 
swaps the states \k) and \l) while leaving all other states 
untouched. Performing now a sequence of mixing oper- 
ations using all possible combinations k < I e {l,...,d} 
for Uki and P = 1/2, one gets rid of the remaining off- 
diagonal elements of the form \'P~j){fJj \ without changing 
the diagonal ones. Thus p can be made diagonal in the 
basis (|3g) by a sequence of bi-local operations. 



B. Mixing of Ha{d) 

Let Hig be the subspace spanned by {|</',i^j)} for fixed 
iQ. In a first step we will show that one can depolarize 
all subspaces Hi„ independently, while in the second step 
we are going to prove that these subspaces can be mixed 
with each other. 

To depolarize Tii, one just has to keep the state 
|1) and performs (randomly) one of the cyclic permu- 
tations of the states {|2), ..., |d)} each with probability 
p = jzr[- Similary, one depolarizes 7^2 by keeping the 
states {|1), 1 2)} and performing with probability p — -j^ 



one of the cyclic permutations of the states {|3), ..., jd)}. 
Since TCi is already depolarized, it is not affected by this 
operation. One can continue in the same way until one 
has depolarized Tid-i- So the antisymmetric part of the 
density operator has now the form 



d-l 



AdpAd^J^'^^ Yl l'^y)('^i; 



(42) 



3=i+l 



The second step starts by mixing of 7id-i with Ti.d-2, 
i.e. equalizing the coefficients ad-i and ad-2- To achieve 
this, both Alice and Bob swap the states \d — 1) and 
|d— 2) with probability's = |i oi' both apply the identity 
operator with probability I ~ ps- If one now depolarizes 
Tid-i and TLd-2 independently as described in step 1, one 
finds that \^-^_^^^){^-^_^^^\ and |'/'(rf_2)j)('^(d-2)jl ^^^*^ 
all the same weight now, i.e. the coefficients are equal. 
Thus the subspace TLd-i is completely mixed up with 
'Hd-2- One now continues by mixing the subspaces TLd-s, 
with {7id-2,Wd-i} and so on, until one reaches Tii. 

We investigate now one particluar 

step in this procedure, namely the mixing of Tifc with 
{Hk+i,'Hk+2, ■■■■I'^d-i}- Both Alice and Bob swap the 
states \k) and |fc+l) with probability's — d'^^k+i "-"^ both 
apply the identity operator with probability 1 — Ps- After 
this one depolarizes TLj (j < k) independently, then mixes 
7id-i with Ti.d~2 as described above. Next one mixes 
Ti-d-d with {Ti.d-2,'Hd-i} and continues in this way until 
one has mixed Hk and {Hk+i,'Hk+2, ■■■■,'Hd-i} ■ It can 
be checked that after this procedure all weight factors dj 
are equal for {j < k). 

So once one has mixed Hi with {7^2,^3, ....,7id-i}, 
one has achieved that the whole antisymmetric subspace 
is completely depolarized, i.e. it can be written as UdAd- 



C. Mixing ofH^d) 



Note that the depolarization of the antisymmetric sub- 
space also mixes the subspace spanned by {|<y3^)} in a 
similar way. Here we show now that one can also depo- 
larize the subspace spanned by {|xfe} and finally that also 
these two subspaces together can be mixed up. To prove 
the first step, let us define the operation {// as follows: 



Ui\k >= \{k + l)mod{d) > 



(43) 



Performing now Ui i^ Ui,l ^ {I, ....,d} each with proba- 
bility p=^ ensures that the subspace spanned by {Ixfc)} 
is completely depolarized, i.e. that p has now the form 

d d 

p = aAd + b J2 H){vtj\+^J2\^k}{Xk\. (44) 

For the second step, we define the unitary operation T: 



T\j) 



1 



k=l 



,i27r(j-l)(fc-l)/d 



\k). 



(45) 



One can check that if we perform the operation T ®T 
with probabihty pt — -ttt and the identity operation 
with probabihty pi — 1 —pT, the diagonal elements of the 
symmetric subspace will each be identical to — did+i) ■ 
The introduced off-diagonal elements can be eliminated 
using the procudure explained above. Note that the an- 
tisymmetric subspace is not affected by this kind of op- 
erations and will thus remain untouched. So finally we 
managed to show that p can be converted to the standard 
form (12) by a sequence of local operations (p7|). The 
weight in Tia (d) was not affected by any of the used mix- 
ing operations, which ensures that ti (Adp) — ti(AdPa)- 



APPENDIX C: BETTER BOUND FOR /3 

In this Appendix we prove that for any vector |^) of 
the form (y|) 



(*|(Q-/3P)®^|*)> 
for /3 < Pn, where 

/3jv 



3^ V^ /3. 



XN 



3JV/37VI/3 ' 



(46) 



(47) 



and xn — a;*(l — 0{1/N)) as A^ — > 00, where x* = 
3(1 - 3"^/3)i/3, and 0(l/iV) denotes a quantity of the 
order of 1/A^. 

The proof is by induction. For A^ = 1, we have /3i = 
1/2. Assuming that the statement holds for A^ — 1, we 
observe that 



1 



vI'|Q(Q-/3P)«^^-^|vI/)>— 1- 



P 



P- 



N-l 



(48) 



and the same holds for all possible permutations with 
respect to the copies. Adding LHS's, dividing the sum 
by A^, and rearranging various terms, we find for large A^ 



(*|(Q-^^/3P)®^|*)> 



\N 



1 



N 

JL 

Pn- 



(49) 



NP^ 
~27" 



il + 0{l/N))-0{p'). 



For P ~ 0{Pn), the last correction term in the above 
inequality can be safely neglected. It is easy to check 
that the right hand side has zero at /3 = Pn, i.e. that 



1 



M,\(Q^PP)^''\^.)>— I 



JL 

Pn 



(50) 



for P < Pm] the statement holds thus for every A^. 

The above result provides a better bound for A^- 
undistillability: the states pa such that 



P< 



3Ar/37vi/3 
are A^-undistillable. 



(1-0(1/^)) 



(51) 
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FIG. 1. Separability and distillability properties of pa (|l2|). 



